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Abstract 

Understanding hadron structure within the framework of QCD is an ex- 
tremely challenging problem. In order to solve it, it is vital that our thinking 
should be guided by the best available insight. Our purpose here is to explain 
the model independent consequences of the approximate chiral symmetry of 
QCD for two famous results concerning the structure of the nucleon. We show 
that both the apparent success of the constituent quark model in reproduc- 
ing the ratio of the proton to neutron magnetic moments and the apparent 
success of the Foldy term in reproducing the observed charge radius of the 
neutron are coincidental. That is, a relatively small change of the current 
quark mass would spoil both results. 
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The chiral properties of QCD have been the subject of considerable attention, from 
chiral quark models [|l]-0] to the less ambitious, but more systematic, approach of chiral 
perturbation theory [|J. Most recently one has begun to realize the importance of chiral 
symmetry in describing the dependence of hadron properties like masses || and magnetic 
moments || on quark mass. This is vital if one is to compare lattice QCD calculations, 
which are presently confined to current quark masses, m, of order 40-80 MeV or higher, 
with experimental data. 

For our purposes the essential point is that chiral symmetry is dynamically broken. The 
resulting Goldstone bosons enter the calculation of hadron properties through loops which 
lead to a characteristic dependence on m which is not analytic. Indeed for the magnetic 

1 /2 

moment of the nucleons one finds a leading non-analytic behavior proportional to m g . In 
the chiral limit oc m q and 

yu p = Hq ~ am n + 0(ml), 

/^S + ^ + oK). (i) 

It is a crucial property of the leading non analytic (LNA) coefficient, a, that it is entirely 
determined by the axial charge of the nucleon and the pion decay constant (both in the 
chiral limit): 

„ 9a m n , . 

Taking the one-loop value of 5u(= Fx + D\ — 0.40 + 0.61) from chiral perturbation theory 
we find a = 4.41. (Note that all magnetic moments will be in nuclear magnetons (jj,n) 
and all masses in GeV). 

Clearly the LNA term is large, of order 0.6 /i^r, at the physical pion mass. This is one 
third of the magnetic moment of the neutron. Provided the G(m%) terms are small at the 
physical pion mass we can use Eq.(|I|) to extract the proton and neutron magnetic moments 
in the chiral limit: 



$ = /i n - am^ s . (3) 

One then finds a model independent expression for the dependence of the proton to neutron 
magnetic moment ratio on the pion mass: 
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Constraining the chiral expansions to reproduce the experimental proton moment yiP and 
the experimental ratio /x p /|/i n | provides 

$ = 3.41 ft*, A = 1.37, (5) 

and 

;/ " 1.37 + 0.09^- + OK). (6) 



As a consequence of Eq.@, we see that the ratio of the p to n magnetic moments varies 
from 1.37 to 1.55 (a variation of order 13%) as m n varies from to 2m^ hys . In terms of the 
underlying quark mass, such a variation corresponds to a current quark mass variation from 
to just 20 MeV. Within the constituent quark model this ratio would remain constant at 
3/2, independent of the change of quark mass. 

A study by Leinweber et al. || suggests a model independent method for describing the 
mass dependence of baryon magnetic moments which satisfies the chiral constraints imposed 
by QCD. We briefly summarize the main results of that analysis. A series expansion of /x p ( n ) 
in powers of m n is not a valid approximation for m n larger than the physical mass. On the 
other hand, the simple Pade approximant 

p(n) 

nPM = ^0 (7) 

Mo 

has the correct leading non-analytic (LNA) behavior of chiral perturbation theory 

and also builds in the expected behavior at large m n . At heavy quark masses we expect 
that the magnetic moment should fall off as the Dirac moment 




for the proton (upper) and neutron (lower) to the chiral limit. The curves are constrained to pass 
through the experimentally measured moments which are indicated by asterisks. 



as to^ becomes moderately large.Q A fit of the Pade approximant to lattice QCD data |||| 
leads to predictions of the magnetic moments of 2.90(20) and —1.79(21) fx N to be compared 
with 2.793 and —1.193 /ijy for p and n respectively. 

Figure [I] shows a similar fit to the lattice data, this time constrained to pass through 
the experimental moments, and providing the solid curve in Fig. ^ for the p/n ratio of 
magnetic moments. The Pade approximate fit parameters are (no, (3) = (3.33,0.527) and 
(—2.41, 0.427), for p and n respectively. Figure |^ also shows the the result of the constituent 
quark model (dashed line) and the variation of the ratio predicted by the leading non-analytic 
behavior of chiral perturbation theory in Eq. (§) (dotted line). The importance of the terms 
of order m^. and higher are revealed by the ratio calculated using the Pade approximant 



1 Note that this form is valid provided cx m, which seems to be true for up to at least 1 
GeV within lattice simulations. 
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FIG. 2. Ratio of the magnitudes of the proton to neutron magnetic moments. The solid curve 
describes the predictions of the Pade approximant while the dashed line denotes the constituent 
quark model prediction of 3/2. The dotted line is the leading non-analytic behavior of chiral 
perturbation theory. The experimental measurement is indicated by the solid point. 



of Eq.fl7|) (solid curve). The values of /Zq var y slightly in the chiral expansion and the 
Pade due to these small higher order corrections at the physical pion mass. However, it is 
important to note that the slopes of the curves agree exactly in the chiral limit, as demanded 
by chiral perturbation theory. 

The key point is that the ratio displays a significant quark mass dependence. It is roughly 
linear in until m n is of order 2mP hys . It is amusing to imagine the excitement had the 
pion mass been 100 MeV heavier at 240 MeV where the Pade crosses the constituent quark 
model prediction of 3/2. However the constituent quark model prediction really corresponds 
to the m n — > oo limit, and Fig. || suggests this limit is approached rather slowly. 

The surprising consequences of chiral symmetry for this famous ratio naturally leads us 
to reconsider the neutron charge radius. The squared charge radius of the neutron (< r 2 >™ h ) 
is obtained from the slope of the neutron electric form factor, Gb„(Q 2 ) as Q 2 — > 0: 
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< r 2 >1= -Q^- 2 G E (Q 2 )\ Q , =0 (8) 

The Sachs electric and magnetic form factors can be written in terms of the covariant vertex 
functions F% and F 2 as 



G E (Q 2 )=F 1 (Q 2 )-^-F 2 (Q 2 



G M {Q 2 )=F 1 {Q 2 ) + F 2 {Q 2 ). (9) 

Note that for a neutral charge particle Fi(Q 2 = 0) vanishes and hence F 2 (Q 2 = 0) is simply 
the magnetic moment of the particle. Now the charge radius squared of the neutron can be 
written as 

< r2 >^- 6 W^" iQ2)l ^ + 2-k- <10) 



Experimentally < r 2 >™ h = -0.113 ± 0.003 ± 0.004 fm 2 0, while the last term in Eq.flK 
the Foldy term [11|, is numerically —0.126 fm 2 . 



The close agreement between the Foldy term and the observed mean square charge radius 
of the neutron has led to considerable controversy. It has been argued that the difference, 
namely the term involving the Dirac form factor (i*i n ), should be interpreted as the true 
indication of the intrinsic charge distribution of the neutron. Clearly this would be quite 
insignificant. On the other hand, decades of modeling the structure of the nucleon have 
suggested that the neutron must have a non-trivial intrinsic charge distribution. Pre-QCD 
it was clear that the long-range tail must be negative, corresponding to the emission of a 
negative pion (n — > p7i~), but old-fashioned meson theory was incapable of describing the 
interior of the neutron. Post-QCD this was resolved in the cloudy bag model PJjj], where 
the convergence of an expansion in numbers of pions was assured — provided the quark 
confinement region was fairly large and the decuplet states (in this case the A(1232)) was 
included on the same footing as the nucleon ||13|| . The neutron charge distribution then 
originated mainly from the Fock component of its wave function consisting of a 7r _ cloud 
and a positive core of confined quarks. Alternatively, within the constituent quark model, it 
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was proposed that the repulsive gluon exchange interaction between the two d-quarks would 
tend to force them to the exterior of the neutron — again yielding a positive core and a 



negative tail fll4 |. 



In view of these expectations of an internal charge distribution, the interpretation of 
< r 2 >™ h in terms of the Foldy term has been controversial. Isgur has recently shown that 
a careful treatment of relativistic corrections for the calculation of < r 2 >™ h , in a quark- 
di-quark model, leads to a recoil contribution that cancels the Foldy term exactly, hence 



restoring the interpretation in terms of an intrinsic charge distribution — see also |T6 
We now show that the study of the chiral behavior of < r 2 >™ h and fi n supports this 
idea, establishing in a model independent way that the observed similarity between the 
experimental value and the Foldy term is purely accidental. 

It is a little appreciated consequence of the approximate chiral symmetry of QCD that 
the mean square charge radius of the nucleon has a leading non-analytic term proportional 



to lnm-Tr HT7 



2 . p(n) I _ _ 1 + % 



< r >2' Ilna = =F 



where the upper and lower sign correspond to p and n respectively. As a result, the charge 
radii of both p and n diverge logarithmically as the quark mass tends to zero. Physically 
this is easy to understand; as — * the Heisenberg Uncertainty Principle allows the pion 
cloud, and therefore the charge density, to extend to infinite distance. For the magnetic 
moment, on the other hand, there is no divergence — indeed the neutron magnetic moment 
increases in magnitude by about 30% as the pion mass moves from its physical value to 
zero. (Loosely speaking, even though the pion may be at a large distance it moves slowly; 
its angular momentum is constrained to one by angular momentum conservation.) 

To summarize, whereas a change of order 5 MeV in the light quark mass leads to a 30% 
change in the Foldy term, the neutron charge radius < r 2 >™ h becomes infinite! Hence, the 
similarity of < r 2 >™ h and the Foldy term is purely an accident. A small change in the quark 
mass leads to completely different values. This physics is not captured in the constituent 



quark model where a 5 MeV change in the light quark mass corresponds to a change in the 
constituent quark mass from roughly 340 to 335 MeV. In this case the neutron charge radius 
originates in the one-gluon-exchange interaction which is proportional to the inverse square 
of the constituent quark mass and therefore < r 2 >™ h would change by only 3%. 

In summary, chiral perturbation theory provides model independent constraints on the 
quark mass dependence of nucleon magnetic moments and charge radii which compel one 
to conclude that the apparent success of the constituent quark model to predict the p/n 
magnetic moment is accidental. Had the pion mass been lighter than the observed value, 
the p/n ratio would drop further from the constituent quark model prediction of 3/2; the 
latter corresponding to the m w — > oo limit. The coincidence of the Foldy term and the 
observed neutron charge radius is also accidental. Here a small change in the quark mass to 
the chiral limit increases the neutron moment by about 30% while the charge radius becomes 
infinite. These results, which are a rigorous consequence of the chiral symmetry of QCD, 
cannot be simulated in conventional constituent quark models. 
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